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Abstract 

Short baseline neutrino experiments such as LSND and MiniBooNE seem to suggest the ex- 
istence of hght sterile neutrinos. Meanwhile, current cosmic microwave background (CMB) and 
big bang nucleosynthesis (BEN) measurements place an upper bound on the effective number of 
light neutrinos, N^jf and the PLANCK satellite will measure N^ff to a much higher accuracy 
and further constrain the number of sterile neutrinos allowed. We demonstrate that if an MeV 
dark matter particle couples more strongly to electrons and/or photons than to neutrinos, then 
p-wave annihilation after neutrino decoupling can reduce the value of N^f / inferred from BEN and 
PLANCK. This mechanism can accommodate two eV sterile neutrinos even if PLANCK observes 
Nf,ff as low as the standard model theoretical value of 3.046, and a large neutrino asymmetry is 
not needed to obtain the correct primordial element abundances. The dark matter annihilation 
also weakens the cosmological upper bounds on the neutrino masses, and we derive a relationship 
between the change in these bounds and the corresponding change in -/Ve//. Dark matter with an 
electric dipole moment or anapole moment is a natural candidate that exhibits the desired prop- 
erties for this mechanism. Coincidentally, a dark matter particle with these properties and lighter 
than 3 MeV is precisely one that can explain the 511 keV gamma-ray line observed by INTEGRAL. 
We show that the addition of two eV sterile neutrinos allows this kind of dark matter to be lighter 
than 3 MeV, which is otherwise ruled out by the CMB bound on N^fj if only active neutrinos are 
considered. 
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I. INTRODUCTION 



The standard three-neutrino paradigm has been very successful in describing the oscil- 
lation phenomenon associated with solar, atmospheric, reactor and accelerator neutrinos. 
However, in recent years, there have been some hints of deviations from this three-neutrino 
paradigm in short baseline neutrino experiments such as LSND [l| and MiniBooNE 

aa. 

In LSND, the transition probability of — t- Ue was measured, but the data indicated 
a 3.8(7 excess of Ue events. Since only antineutrinos were involved in LSND, it was then 
noticed that if CPT is violated in the neutrino sector such that the mass-squared splittings 
for neutrinos and antineutrinos are different, all the data could be reconciled 41. But this 

ri 

scheme was contradicted by the KamLAND [51 data. Another CPT- violating scheme was 

a, 

neutrmo data |3 . 

n 

Meanwhile, the idea of sterile neutrinos [g] has also been invoked to resolve the LSND 
anomaly. But it turns out that adding only one eV sterile neutrino to the standard three- 
neutrino picture, namely the (3+1) model, is not enough The (3+2) model (with two 
eV sterile neutrinos) was proposed and shown to fit the data much better than the (3+1) 



model 



io|. 



The purpose of MiniBooNE was to confirm or exclude LSND. In the first release of data, 
they found no evidence for an excess of in the u^j — )■ Uf, search [2:]. But their later search 
for z/^ — )■ Ue appeared to be consistent with LSND jsj]. In contrast to their previous z/^ — > z/g 
search, MiniBooNE's latest result indicates that the neutrino and antineutrino data are 
consistent with each other 



(such as CDHS 



11| . However, the tension between appearance and disappearance 



12| and Bugey 13[) short baseline experiments remains l4|. A (3+3) model 



has been considered by pJi] , but the global fit to appearance and disappearance experiments 
is still poor. 

Recently, a further hint of sterile neutrinos has emerged from a theoretical re-evaluation 

n 

of the expected mean reactor antineutrino flux emitted from nuclear reactors UM. The new 

n 

prediction suggests a flux 3% higher than what was previously assumed [17[. This implies 
that all reactor neutrino experiments with a baseline smaller than 100 m have measured a 
deficit in the Ue events. Motivated by this new observation, the global short baseline neutrino 
oscillation data has been re-analyzed in the (3+1) and (3+2) models 18|. It is shown that 
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while the (3+1) model is still insufficient, the global fit improves significantly in the (3+2) 
model. The conclusion is that the short baseline neutrino experiments now favor two sterile 
neutrinos with the best-fit masses being 7714 = 0.68 eV and ms = 0.94 eV 18|. 

While short baseline neutrino experiments suggest the existence of eV sterile neutrinos, 
there are useful constraints from cosmology. For instance, the bound from Big Bang Nucle- 
osynthesis (BBN) is quite stringent. The existence of extra relativistic degrees of freedom, 
if thermalized, will increase the Hubble expansion rate and cause the weak interactions to 
freeze out at earlier time (higher temperature). This will increase the neutron-to-proton 
ratio and may result in an overproduction of primordial helium and deuterium. The recent 
analyses lead to N^ff < 4.26 (loj and N^ff < 4.1 |2oj at 95% C.L. limit. The best fit derived 
by 19| is Neff = 3.86. As a result, while one fully thermalized eV sterile neutrino is pre- 



ferred by BBN, two fully thermalized eV sterile neutrinos are disfavored. This bound could 
be circumvented if there is a positive z/g degeneracy parameter (or equivalently, neutrino 
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22 1 . With a non-zero Up de- 



asymmetry), ^e, which reduces the neutron-to-proton ratio 
generacy parameter, the neutron-to-proton ratio is given by n/p = exp{ — (m^ — mp)/T — ,^e}. 
So the purpose of the degeneracy parameter is to reduce n/p which has been elevated by 
the two thermalized eV sterile neutrinos. In order to accommodate two eV sterile neutrinos, 

one would need the degeneracy parameter to satisfy 0.03 < < 0.14 at 95% C.L. with 

! I 

the best-fit value being = 0.064 [19]. 

In fact, if is really as large as approximately 0.1, the active-sterile mixing will be 
delayed through the suppression of the medium mixing angle j2]| . In this case, the medium 



mixing angle wi 
asymmetries 



1 be suppressed by the matter effect which is dominated by the neutrino 
23j |. The active-sterile transition probability is even suppressed by the so- 
called quantum Zeno effect near the MSW resonance |2^]. Thus, the resulting effect is that 
the active-sterile mixing will occur only after neutrino decoupling and the sterile neutrinos 
will not be thermalized during the epoch of BBN. This understanding has recently been 
confirmed by {25 1 which solves the full quantum kinetic equations 26 1. As a consequence, 
the effective number of neutrinos remains at N^ff = 3.046. Without the thermalized sterile 
neutrinos, the existence of would reduce n/p and the question is how much can be 
tolerated by BBN. It turns out that the answer depends on the value of the mixing angle ^13. 



Recently, DAYA-BAY 



27j and RENO |28j have observed sin^ ( 2^13 ) = 0.092+0.016 (stat.) + 



0.005 (syst.) and sin^ ( 2^13 ) = 0.113 + 0.013 (stat.) + 0.019 (syst.) at 68% C.L. respectively. 
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Taking into account all of the measured mixing angles for the three active neutrinos and 
;he effect of neutrino oscillations, it has been shown that ~ 0.1 is indeed allowed by BBN 



29|. 



The CMB observations 
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j32| also place bounds on the effective number of neutrino 
degrees of freedom N^ff. (See [33| for a discussion of the effect of N^ff on the CMB fluctu- 
ations.) At recombination, the temperature of photons and hence their energy density is 
extremely well-measured. This means that the cosmic neutrino background with tempera- 
ture Ty = (4/ll)^/^Ty is well-predicted accordingly. Assuming slight heating of the neutrinos 
due to e^e~ annihilation, the theoretical prediction for the effective number of neutrinos is 
N^ff = 3.046 
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35|. Although the neutrino energy density cannot be measured directly. 



it can be deduced from the CMB measurements. The knowledge of neutrino energy density 
leads to the bound on iVe//. 

The values of N^ff from a combined analysis of WMAP, baryon acoustic oscillation 
(BAG) and Hubble constant (Ho) yields N^ff = 4.34+°;^^ (68% CL) [sO]. The further 
addition of the Atacama Cosmology Telescope [31[ and the South Pole Telescope [32| lead 
to Neff = 4.56 ± 0.75 (68% CL) and N^jf = 3.86^0.42 (68% CL) respectively. Combined 



datasets have been used by Archidiacono et al. [361 derive N^ff = 4.08_o:68 (95% CL). 
However, these constraints on N^ff are derived under the assumption that the additional 
particles are massless. Taking into account of the masses of two light sterile neutrinos, it has 
been shown that ^nii, < 0.9 eV 37L which is clearly in conflict with the best-fit masses 



m4 



0.68 eV and mg = 0.94 eV 



18| derived from neutrino experiments. Nevertheless, 



this issue appears far from settled. For instance, in a recent analysis, it has been shown 
that these best-fit masses are marginally consistent with the cosmological data 38|. But 



two other recent analyses conducted by 



39[ and j40| seem to suggest ^m^j < 0.6 eV (with 



N^ff = 3.58) and Xl'^i^ ~ 0.556 eV (with N^fj = 3.839) respectively, while Ref. j41|] obtains 
S rn^< 0.32 eV and Ref. j42| gives J2 "^^^ ^ 0.34 eV. On the other hand, the remarks made 
by |43| are illuminating — they show that whether the cosmology of the (3+2) model is 
viable depends significantly on the choice of datasets included and the ability to control the 
corresponding systematic uncertainties. They thus conclude that it is premature to claim 
that two eV sterile neutrinos are either cosmologically favored or ruled out. As a result, it 
appears that we need a better understanding of the systematic uncertainties in the analyses 
of cosmological datasets in order to nail down an indisputable bound on '^m^. 
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To recapitulate, we have just reviewed the current status of sterile neutrinos in the context 
of both particle physics and cosmology. The forthcoming data from the PLANCK satellite 
will improve the bound on iVg// to an accuracy of 0.20 at 2(j. For instance, if PLANCK 
sees Nf.ff < 4.60, then it is unclear how to reconcile the CMB bound with two eV sterile 
neutrinos. In a worse case, if PLANCK sees A^e// ^ 4.0, it would appear that two eV sterile 
neutrinos are completely ruled out. One of the main purposes of this paper is to propose 
a novel mechanism, which involves MeV dark matter p-wave annihilation, to accommodate 
two eV sterile neutrinos despite the stringent CMB bound. We will demonstrate that this 
mechanism is capable of accommodating two eV sterile neutrinos even if PLANCK observes 
Nf.ff as low as the theoretical value of 3.046. It can also accommodate the BEN bound on 
the number of sterile neutrinos without introducing a large neutrino asymmetry. We will 
describe this mechanism in Sec. |TT] and Sec. IIIII Then we will identify a few natural dark 
matter candidates that exhibit the required properties for this mechanism in Sec. IIVI 

Interestingly, a dark matter particle with precisely these required properties and lighter 
than 3 MeV can explain the 511 keV gamma-ray line observed by INTEGRAL 4^. In Sec. 
IVl we will show that two eV sterile neutrinos allow this kind of dark matter to be lighter 
than 3 MeV, which is otherwise ruled out by the CMB bound on A"ejj if there are no sterile 
neutrinos. 



II. REDUCING iVe// THROUGH MEV DARK MATTER ANNIHILATION 

In this section, we present a novel mechanism to reduce Nf,ff through MeV dark matter 
annihilation. Assuming only the three active neutrinos, this mechanism has recently been 
invoked by j45| to place a lower bound on the mass of MeV-scale dark matter. As we will 
see, if eV sterile neutrinos exist, the bound derived in {45 1 will be somewhat relaxed. 

The main idea behind this mechanism is the following. An MeV dark matter particle 
that couples more strongly to electrons and/or photons than to neutrinos will heat up the 
electron-photon plasma through pair-annihilation when it becomes non-relativistic before 
its abundance freezes out. If this occurs after neutrino decoupling, then the ratio of the 
neutrino temperature to the photon temperature will be reduced, leading to a decrease in 
Neff- This process is similar to the heating that results from the annihilation of electrons 
and positrons when they become non-relativistic. 
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The fact that a generic hght particle, which remains in thermal equilibrium with the 
photons throughout the epoch of BBN, will reduce iVg// was first briefiy discussed by 46 1 
in the context of BBN. They argued that this effect may lead to an underproduction of 
primordial helium and deuterium. More recently, the effect of light dark matter annihilation 
on BBN has been investigated by [47i] . 

MeV dark matter has been invoked to explain the observed 511 keV 7-rays by INTEGRAL 



48l | and the cosmic 7-ray background at 1 — 20 MeV 49|. It has also been found to have 



interesting effects on large-scale structure jsO]- A supersymmetric model with the lightest 
supersymmetric particle of MeV scale has been proposed in [51]. Then, in the context of 
supersymmetric models with gauge mediated SUSY breaking, it has been shown that MeV 
dark matter serves as an example of WIMPless candidate that naturally leads to the correct 
relic abundance [52I. 

A constraint on MeV dark matter from core-collapsed supernovae has been derived by 



531] . But their bound can be evaded if the scattering cross section between the dark matter 



and neutrinos is negligible. Since we are interested in MeV dark matter that couples more 
strongly to electrons and/or photons than to neutrinos, this bound is not relevant to our 
consideration. 

Meanwhile, dark matter annihilation (which generates the correct relic abundance) near 
the epoch of recombination will cause distortions to the CMB fluctuation spectrum and so a 
lower bound on the dark matter mass can be derived |5J] . (Notice that dark matter annihi- 
lation also distorts the CMB spectrum jssl], although the bound is much weaker.) This effect 
excludes dark matter with masses < 1 — 10 GeV, which is apparently in conflict with MeV 
dark matter. However, this bound is only applicable to s-wave annihilation, for which (erf) 
remains constant between the dark matter freeze-out and the epoch of recombination. For 
j9-wave annihilation which is velocity-dependent, the annihilation rate at the recombination 
epoch is negligible, and so this bound can be evaded. Therefore, we require an MeV dark 
matter particle with p-wave annihilation for our mechanism of reducing N^ff. 

The following discussion is based on Ref. 45|]. Let XX be the pair of dark matter 
particle and anti-particle. (Notice that for self-conjugate and non-self-conjugate scalars, the 
notation xX really means xX ^-^id xX* respectively.) Consider the case in which the dark 
matter annihilates entirely after the neutrinos have decoupled in the early universe. This 



occurs at a temperature of ~ 2 — 3 MeV [3J, l56| . The extent of the heating due to dark 
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matter annihilation can be derived from entropy conservation 57|, |58| . The total entropy 
before XX annihilation is proportional to 



S = ^ { Pe+e- + P7 + PXX + Pe+e- + Pj + PxX 



(2.1) 



while after xX annihilation, it becomes 

S =^ (Pe- 



+ p-y+ Pe+e~ +P-y) ■ 



(2.2) 



For a relativistic particle, the pressure is related to the density by p = p/3, so the total 
entropy density can be written as 58l |: 

Ptot + Ptot 2 vr^ 



(2.3) 



T 45 

where g^:S is the sum of total bosonic spin degrees of freedom and 7/8 times the total 
fermionic spin degrees of freedom. Then, the total entropy equals to 

2 7r2 



S 



45 



■9*s {RTY 



(2.4) 



which is conserved during the process where any particle species becomes non-relativistic 
and pair-annihilates. Thus, the ratio of the values of RT before and after xX annihilation 
is given by 



{RT), 



9*Sf 



1/3 



(2.5) 



where (7^,5^ and g^Sf are respectively the values of g^s for the relativistic degrees of freedom 
in thermal equilibrium before and after xX annihilation. 

If the XX annihilation occurs after neutrino decoupling, then the neutrinos do not share 
in the heating from the annihilation. This means that RTy remains constant before and 
after xX annihilation and T;^ oc R~^ . But the photons and electron-positron pairs are heated 
up according to Eq. (12.51) . Therefore, the ratio of T^, to Ty after xX annihilation will be: 

f ^ .1/3 



(7/8)4 + 2 



XX 



(7/8) 4 + 2 + (7/8) 9 
1/3 



(7/8)4 + 2 



for a fermionic x '■, 



for a bosonic x ■ 



(2.6) 



_(7/8)4 + 2 + g_ 

where g is the total internal degrees of freedom for the xX pair. For a self-conjugate scalar 
dark matter particle, we have g = 1, while for a non-self-conjugate scalar dark matter 



particle, we have g — 2. Further, for a spm-1/2 Majorana dark matter particle which is 
self-conjugate, we have g — 2; while for a spin-1/2 Dirac dark matter particle, we have 
g — 4:. Subsequent e'^e~ annihilation further heats up the photons, so the resultant ratio of 

Ti, to Try is given by 

T / zL \ / T \ 

(2.7) 



In terms of iVe//, the energy density of all relativistic neutrinos can be written as 

The CMB observations constrain p^j at fixed and so a change in T^^jT^ with respect to 
its conventional value (4/11)^/^ (due to e"'"e^ annihilation) would be interpreted as a change 
in Ai'e//. But as we have shown above, XX annihilation reduces the value of T^/T^ by a 
factor of {Ti,/T^)^^ relative to (4/11)^/^. So the apparent value of N^ff inferred from CMB 
observations will be a factor of (Tj^/T^)^^ smaller than it would without xX annihilation. 

In this calculation, we have neglected the effect of the annihilating dark matter particle on 
the expansion rate when the neutrinos drop out of thermal equilibrium (at which point the 
dark matter is relativistic), as well as the possible effect of two additional sterile neutrinos. 
Both of these will increase the expansion rate and cause the neutrino decoupling temperature 
to increase. We can estimate the resultant effect as follows. The Hubble rate is proportional 
to ^/g^ T"^, while the weak interaction rate is proportional to T^. The neutrino decoupling 
temperature 7^ is defined as the temperature at which the weak interaction rate drops to 
roughly the same magnitude as the Hubble rate. This implies that oc ( g*s)^^^ and so the 
modified neutrino decoupling temperature will be given by 

^modified _ / 2 + (7/8) 4 + (3 + 2) (7/8) 2 + (7/8) g \ 

^original 3 + (7/8) 4 + 3 (7/8) 2 y ' ^ ' ' 

where, by writing the factor {7/8)g in the above expression, we have assumed a fermionic 
dark matter. But similar arguments apply to a bosonic dark matter. As an example, one 
can take g = 2 for a Majorana dark matter. Then, we will have Tyodificd/jiorigmai _ qj 
which indicates that the change is small. Thus, we will ignore this effect here and in the 
following section. 
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III. EV STERILE NEUTRINOS AND MEV DARK MATTER 



What we have considered in the last section was the effect of XX annihilation occurring 
completely after neutrino decoupling. In reality, it is possible that xX become non- 
relativistic and start to annihilate before neutrino decoupling. In this case, the neutrinos 
and photons will share in the heating from xX annihilation before neutrino decoupling, while 
the residual xX annihilations after neutrino decoupling will heat the photons alone. 

Our treatment follows Ref. 45|]. After neutrino decoupling, the xX annihilation continues 
to heat up the photons relative to neutrinos until the xX particles drop out of thermal 
equilibrium. But the exercise in the last section reveals that the ratio of to in each 
step of XX annihilation depends solely on g^Sr and g^Sf To quantify this continuous effect. 



we would need to consider the quantity I{T^) given by (see [57| for a similar calculation): 
S 1 

_nn\g r . f / 2 , r ^rr. , 

(3.1) 

where the integration variable is x = p-^/T^. The "+" sign and "-" sign in the integrand 
correspond to fermionic and bosonic dark matter respectively. Up to the constant factor 
2 7r^/45, the physical meaning of I{T^) is the sum of total bosonic spin degrees of freedom 
and 7/8 times the total fermionic spin degrees of freedom at a given photon temperature 
T^. Notice that in the limit where all particles are relativistic, I reduces to (2 7r^/45) (7*5 
and the integral in Eq. (13. ip simply counts the contribution from xX when they become 
non-relativistic. 

As mentioned above, the xX annihilation will continue to heat up photons relative to 
neutrinos after neutrino decoupling until the xX particles drop out of equilibrium. Thus, 

the ratio Ty to due to accumulative xX annihilation is 

1/3 



(3.2) 



_ (m)\ 

where Tf is freeze-out temperature of the xX particles and is given by Tf ~ as a 



generic result [58[. Since most of the entropy due to xX annihilation is deposited to the 
thermal background around T ~ m^/3 ^ Tf, it is evident from Eq. (13. ip that we can set 
= in Eq. (13. 2p as a good approximation. 



A^eff 



5.0 r 




FIG. 1: Neff as a function of m^/Td for two sterile neutrinos. The curves with g = \, g = 2 and 
g = A correspond to a self- conjugate scalar, Majorana and Dirac dark matter respectively. The two 
dashed lines correspond Nf,ff = 3.046 and N^jf = 4.00. 

We expect each of the two eV sterile neutrinos to contribute a unity to the effective 
number of neutrinos. But as a consequence of XX annihilation, Nejf deduced from CMB 
observation would be given by 

The value of N^fj against m^/Ta is shown in Fig. [H for a self-conjugate scalar {g = 1), 
a spin-1/2 Majorana {g = 2) and a spin-1/2 Dirac {g = 4) dark matter. The curve for a 
non-self-conjugate scalar {g = 2) dark matter appears very similar to that for a spin-1/2 
Majorana dark matter and so it is not shown. 

Note that the small increase in the neutrino decoupling temperature due to the xX P^^- 
ticles and the two eV sterile neutrinos, which we have chosen to neglect, will only serve to 
enhance the heating of the photons relative to the neutrinos, further decreasing N^ff. 

The PLANCK satellite, whose first data will be released very soon, is expected to improve 
the bound on A^e// to an accuracy of 0.20 at 2a. If PLANCK sees A"e// > 5.0, then our 
mechanism of reducing N^ff through dark matter annihilation is not needed to accommodate 
two eV sterile neutrinos. If PLANCK sees 4.0 < Ngff < 5.0, one can easily read off the 
corresponding mass of a self-conjugate scalar, Majorana or Dirac dark matter required to 
accommodate two eV sterile neutrinos from Fig. [1] If PLANCK sees 3.5 < ^^eff ^ 4.0, both 
Majorana and Dirac dark matter can help to accommodate two eV sterile neutrinos, while 
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the self-conjugate scalar dark matter fails to do the job. If PLANCK sees 3.046 < Nefj ^ 
3.5, only a Dirac dark matter can do the job. In fact, a Dirac dark matter is capable of 
accommodating two eV sterile neutrinos even if PLANCK sees A"e// as low as the theoretical 
value of 3.046. 

In addition, we remark that in the unlikely event that PLANCK sees Nf,ff < 3.046 (in 
contradiction to the standard cosmological model), then one possible resolution would be 
that MeV Dirac dark matter annihilation has indeed occurred to reduce A^e//- On the 
other hand, if PLANCK sees the highly unlikely value N^fj < 2.5, then even a Dirac dark 
matter with vanishing mass cannot help to accommodate two eV sterile neutrinos. If this 
happens, two eV sterile neutrinos would be strongly disfavored. One may then conclude that 
either sterile neutrinos are not relevant for the anomalies from the short baseline neutrino 
experiments or some exotic physics needs to come in to rescue. 

As noted above, most of the entropy due to XX annihilation is deposited to the thermal 
background around T ~ m^/3. For an MeV dark matter particle, this means that most 
of the entropy from xX annihilation would have already been released before BBN starts. 
Thus, the value of N^jf inferred from BBN would be roughly equal to that obtained from Eq. 
(13. 3p . In other words, MeV dark matter p-wave annihilation can also help to accommodate 
the BBN bound on the number of sterile neutrinos. This mechanism does not require a large 
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22|. 



neutrino asymmetry, as opposed to the work by 

While we have shown that the annihilation of an MeV dark matter particle can reduce 
Neff to an acceptable value even with two additional sterile neutrinos, we have not yet 
considered the observational limits on discussed in Sec. I. Limits on the neutrino masses 
are, to a first approximation, limits on the total non-relativistic neutrino density, which is 
given by 

py = ^nym^, (3.4) 

where the sum is taken over all (standard model and sterile) non-relativistic neutrinos, 
and ny is the number density of each type of neutrino. If all neutrinos have a common 
temperature, Ty^ then ny is the same for every neutrino and is given by 
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where ( is the Riemann zeta function and C(3) = 1.2. Thus, we can rewrite Eq. (13. 4 p as 

p-(|)2(^)c(3)(|)V=E-- P«) 

If we have an observational bound on p^, at a given Ty, expressed as an upper bound on 
^ m,y, it is clear from Eq. (13. 6 p that a decrease in T^/T^ can compensate an increase in 
^ m,^ and continue to satisfy the observational bound. This implies that the apparent value 
of Y2mi, inferred from cosmological observations will be a factor of (T^/T^)^^ smaller than 
it would without the MeV dark matter annihilation: 



- m ■ p-^) 



cosmology V ' ^ / true \ , 

In other words, the upper bound on ^ ra^ is increased by a factor of {T^jTy)'^^ = I{Td) / /(O). 
It is easy to express this factor in terms of the change in iVe// produced by the annihilating 
MeV dark matter. If N^fj is the value deduced from CMB observations, and N^, is the true 
number of neutrinos (e.g., Ny = 5.046 in the 3 + 2 model), then we see that 

/ true \ / cosmology yiVgjj y 

So the upper bound on J^^i' increased by a factor of (Ny/Neff)^^'^. For the case of two 
sterile neutrinos, Ny = 5.046, and value of A^e// (for a given m^/Td) can be read off of Fig. 
1. Thus, we obtain a relationship between the dark matter annihilation needed to satisfy 
upper bounds on J2 "^^j the prediction for the observed value of N^ff. As an example, if 



the neutrino experiments give "^rriy ^ 1.6 eV [18|, while cosmological observations require 
^ < 0.9 eV [371], then we have (5.046 /A^e//)^/^ > 1-6/0.9, or N^jf < 2.3, which is difficult 
to reconcile with current observations. However, as noted in Sec. I, the situation with regard 
to upper bounds on ^ m^, remains somewhat murky. 

Our simple treatment of limits on ^ rriu as corresponding to an upper bound on p^, alone 
is not entirely accurate. The values of rrii, also determine the temperature at which each 
neutrino becomes non-relativistic, and thereby alter the process of neutrino free-streaming, 
which affects large-scale structure. In our scenario, this process is further changed because 
Ti^/Ty is reduced below its value in the standard model, so that a neutrino of a given mass 
will become non-relativistic earlier than it would without the MeV dark matter annihila- 
tion. However, these effects are below the threshold for detection by current cosmological 



observations and would be difficult to detect even in future surveys 
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Clearly, our scenario differs from the standard cosmological model in predicting a different 
relationship between Pu{T.y) when the neutrinos are highly relativistic and Pu{T^) when they 
are non-relativistic, since the former scales as {Ti^/T^Y while the latter scales as {T^/T^Y. 
This predicted difference could be detectable by future cosmological observations. 



IV. RELEVANT DARK MATTER MODELS 



Since we are interested in dark matter candidates that couple more strongly to electrons 
and/or photons than to neutrinos, the prototypical example could be a Dirac dark matter 
that interacts with the ordinary matter through an electric dipole moment (EDM) or a 
magnetic dipole moment (MDM) 60|: 

C = x{i^-m^u)x + ^^X^^''F,,X + i^X^^''F,,X, (4.1) 

where F^^ and F^^ are the electromagnetic field and dual field strengths respectively. In 
order to be consistent with direct detection experiments, it has been pointed out that thermal 
dark matter (with the correct relic abundance) acquiring EDM or MDM must have a mass 
less than 1 - 10 GeV [61 1. 



For dark matter MDM, the annihilation cross sections are 



6l|: 



MDM (^AfCt) TT 



DM 



DM 



which are s-wave annihilations. But as mentioned earlier, s-wave dark matter annihilation, 
which leads to the correct relic abundance, distorts the CMB fluctuation spectrum and the 
corresponding dark matter masses < 1 — 10 GeV are excluded. So MeV dark matter with 
an MDM is not allowed. 



For dark matter EDM, the annihilation cross sections are 6l| : 



^EDM „ {OEofTl 2 

^XX^e+e-^rel = 40^2 ^rel , (4-4) 

DM 



EDM „, (fl'l;'^)^^ f. r\ 

^xx^77^-i = g4^2 • (4-5) 



DM 



While (T™^^^Wrei is s-wave, it is proportional to In the early universe, Wrei is not sup- 
pressed. So it is possible that cr'^?^„+ ^-frei, which is only proportional dominates over 



13 



(y^^^-yV^ei and is responsible for the correct relic abundance. In this case, the bound from 
CMB fluctuation spectrum does not apply. So an MeV dark matter with an EDM is a viable 
candidate for our mechanism of reducing N^fj. 

In addition, a fermionic (Dirac or Majorana) dark matter may also acquire an anapole 



moment (AM) through the following interaction 62l |: 



/:AM = |^X7'^7'x5"i^M- (4-6) 

where is the coupling constant and A is the cut-off scale. This interaction operator 
violates C and P individually, but preserves T. Interestingly, for a Majorana fermion, both 
of EDM and MDM are forbidden, and the only allowed electromagnetic form factor is the 



anapole moment 63|. At the tree- level, the annihilation cross sections are given by [621 ]: 



Y^e+ e- ^rel — TTT^ ^rel ) l^- ' J 



2 (yf^ a 



which is purely p-wave. In contrast to EDM and MDM which interact with external electro- 
magnetic fields, the anapole moment interacts only with external electromagnetic currents 
= d'^F^ij. But the on-shell external photons do not lead to an electromagnetic current. 
This explains why we have cr^|L.^^frei = at the tree-level. It has been shown that dark 
matter with an anapole moment can generate the correct relic abundance and be consistent 
with the direct detection experiments simultaneously for the mass range from MeV up to 



100 GeV j62|. Obviously, an MeV dark matter with an anapole moment is also a viable 
candidate for our mechanism of reducing N^ff. 

V. 511 KEY GAMMA-RAY LINE 

The INTEGRAL satellite has observed a 511 keV gamma-ray emission line from the 



galactic center 44|. This emission line is expected to arise from non-relativistic positrons 
annihilating with the electrons at rest in the galactic bulge. So it is important to understand 
the source of these positrons. 

One natural way to produce these positrons is through dark matter annihilation. For dark 
matter with mass larger than GeV, the injection energies for the positrons would be too high 
to confine to them in the galactic bulge. Also, a substantial fraction of the positrons might 
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have already annihilated before reaching the required non-relativistic energies. The idea of 
1-100 MeV dark matter was then proposed by [48] to explain this 511 keV gamma-ray line. 
Later, it was shown by [64] that the dark matter annihilation process XX ~^ e^e~ is neces- 
sarily accompanied by the electromagnetic radiative corrections (internal bremsstrahlung) 
which lead to the process xX e+e^7. What they have found was that unless < 20 
MeV, the real gamma-rays emitted from the electromagnetic radiative corrections would 
violate the constraints on the gamma-ray flux from COMPTEL and EGRET. By compar- 
ing the gamma-ray spectrum generated from positron inflight-annihilation and the observed 
diffuse galactic gamma-ray data, the positron injection energies were further constrained to 



be < 3 MeV 



65| . Recently, it has been shown that the MeV dark matter explanation is 



consistent with dark matter halo profiles predicted by numerical many-body simulations for 
a Milky Way-like galaxy 66 1. 



For an MeV dark matter, the only kinematically allowed annihilation modes are e"'"e~, 
photons and neutrinos. But to avoid the direct gamma-ray constraint such as the cosmic 
gamma-ray background jfi?! and to maintain the required positron production rate, the 
annihilations into photons and neutrinos are postulated to be suppressed. Therefore, the 
current consensus is that in order to simultaneously explain this 511 keV gamma-ray line 
and produce the correct relic abundance, we need a dark matter with mass < 3 MeV 
which annihilates primarily into e^e~ through p-wave 



The MeV dark matter required could be a scalar or fermion which annihilates primarily 
into e^e~ through a new light gauge boson 48|, |69|. It could also be a dark matter with an 
EDM or anapole moment as mentioned in Sec. IIVI In other words, the dark matter, which 
can help to accommodate two eV sterile neutrinos and is lighter than 3 MeV, is precisely 
the one required to explain the 511 keV gamma-ray line observed by INTEGRAL. 

One the other hand, as shown in 45|, for any dark matter that couples more strongly 



to electrons and/or photons than to neutrinos and annihilates through p-wave, the lower 
mass bound set by Nejf inferred from CMB is at least > 3 MeV. This mass bound 
applies to any dark matter that may be a self-conjugate scalar, non-self-conjugate scalar, 
Majorana fermion or Dirac fermion. So the possibility of explaining the 511 keV gamma-ray 



^ See, however, [68j for possible astrophysical arguments against the dark matter interpretation of the 511 
keV gamma-ray line. 
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line by dark matter with m y < 3 MeV is ruled out by the lower mass bound derived in 45 1. 
However, the derivation in [45| assumed no sterile neutrinos. With two eV sterile neutrinos, 
it is evident from Fig. [T]that the lower mass bound can easily be relaxed. This reopens the 
window for < 3 MeV. 

From Fig. [T], if we choose a reasonable lower bound on N^fj to be N^fj > 3.046, then 
there is no lower mass bound for a self-conjugate/non-self-conjugate scalar or Majorana 
dark matter, while a Dirac dark matter must have > 4 MeV. So a Dirac dark matter, 
including the one with EDM or anapole moment, cannot be a candidate to explain the 
511 keV gamma-ray line. But a self-conjugate scalar or Majorana dark matter may still 
be viable depending on the value of N^jf to be measured by PLANCK. For instance, if 
PLANCK measures 4.0 < ^e// ^ 4.2, only the self-conjugate scalar is a viable candidate. If 
PLANCK measures 3.5 < A"e// < 3.7, only the Majorana dark matter is a viable candidate. 
Otherwise, it becomes impossible to explain the 511 keV gamma-ray line by dark matter 
annihilation. 

Finally, we would like to compare the models that consist of a Majorana dark matter 
particle for explaining the 511 keV gamma-ray line. So far, all the proposed models require 
the existence of a new light gauge boson and need to postulate that the dark matter an- 



nihilations into photons and neutrinos are somehow suppressed |48l . |5l| • This ensures that 
the Majorana dark matter annihilates primarily into e~^e~ . On the other hand, a Majorana 
dark matter with an anapole moment exhibits all of these required properties in a simpler 
and more natural manner j62|. Firstly, Eq. ( 14. 6 p dictates that the dark matter anapole 
moment couples to an external electromagnetic current. The Majorana dark matter annihi- 
lates through the Standard Model photons. So no new gauge boson is required. Secondly, 
since the photons do not couple to the neutrinos, the dark matter annihilation into neutrinos 
is naturally suppressed. Thirdly, according to Eq. (14. Sp . we have cr^|L.^^ Wrei = at the 
tree-level, and so the dark matter annihilation into photons is naturally suppressed as well. 
The dark matter annihilation into e'^e" is the predominant mode and is p-wave. Therefore, 
an MeV Majorana dark matter with an anapole moment is the simplest and most natural 
candidate to explain the 511 keV gamma-ray line. 
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VI. CONCLUSIONS 



In this article, we presented a mechanism for reducing A^e// through MeV dark matter 
p-wave annihilation. This occurs if the MeV dark matter couples more strongly to electrons 
and/or photons than to neutrinos. Consequently, we showed that this mechanism can help 
to accommodate two eV sterile neutrinos even if PLANCK observes A^e// as low as the 
theoretical value of 3.046. It can also help to accommodate the BBN bound on the number 
of sterile neutrinos without introducing a large neutrino asymmetry. This mechanism also 
weakens the upper bound on the neutrino masses. Further, our model (or indeed any model 
which changes the ratio between T^, and T^) produces a different relationship between Pi,{T^) 
when the neutrinos are highly relativistic and Pu{T^) when they are non-relativistic, since 
the former scales as (T^/T^Y while the latter scales as (T^/Tj)^. This modification to the 
standard cosmological model could become apparent (or ruled out) as the observational data 
improve. 

We pointed out that a natural dark matter candidate which exhibits the desirable prop- 
erties for this mechanism could be the one with an EDM or anapole moment. Interestingly, 
in order to explain the 511 keV gamma-ray line observed by INTEGRAL, we need precisely 
a dark matter that enables this mechanism and is lighter than 3 MeV. We showed that while 
this kind of dark matter lighter than 3 MeV is ruled out by the CMB bound on iVg// if only 
active neutrinos are considered, two eV sterile neutrinos reopen the window for < 3 
MeV. Finally, we argued that an MeV Major ana dark matter with an anapole moment is 
the simplest and most natural candidate to explain the 511 keV gamma-ray line, and at 
the same time, can help to accommodate two eV sterile neutrinos through its purely p-wave 
annihilation into e'^e" . 

It is remarkable that the solution for the anomalies in LSND and MiniBooNE (two eV 
sterile neutrinos) is somehow related to the solution for the 511 kcV gamma-ray observed 
by INTEGRAL (MeV dark matter). On the one hand, we need MeV dark matter p-wave 
annihilation to reduce N^ff so as to accommodate two eV sterile neutrinos despite the BBN 
and CMB bounds. On the other hand, two eV sterile neutrinos are useful in allowing a dark 
matter of exactly the same kind to be fighter than 3 MeV. We thus conclude that eV sterile 
neutrinos and MeV dark matter are complementary to each other. 
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